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ABSTRACT 

We present results from Submillimeter Array (SMA) 860-/im sub-arcsec astrometry and multiwave- 
length observations of the brightest millimeter (Si.i mm = 8.4 mJy) source, SSA22-AzTECl, found near 
the core of the SSA22 protocluster that is traced by Lya emitting galaxies at z = 3.09. We identify a 
860- /im counterpart with a flux density of Sseoum — 12.2 ± 2.3 mJy and absolute positional accuracy 
that is better than 0'.'3. At the SMA position, we find radio to mid-infrared counterparts, whilst no ob- 
ject is found in Subaru optical and near-infrared deep images at wavelengths < 1 /jm (J > 25.4 in AB, 
2a). The photometric redshift estimate, using flux densities at > 24 /im, indicates z p hot = 3. 19^q 35, 
consistent with the protocluster redshift. We then model the near-to-mid-infrared spectral energy 
distribution (SED) of SSA22-AzTECl, and find that the SED modeling requires a large extinction 
(Ay ps 3.4 mag) of starlight from a stellar component with M sta r ~ 1O 1O ' 9 M0, assuming z = 3.1. 
Additionally, we find a significant X-ray counterpart with a very hard spectrum (r e g = — 0.34^ 06 y, 
strongly suggesting that SSA22-AzTECl harbors a luminous AGN (Lx ~ 3 x 10 44 ergs s~ 4 ) behind a 
large hydrogen column (TVh ~ 10 24 cm -2 ). The AGN, however, is responsible for only ~ 10% of the 
bolometric luminosity of the host galaxy, and therefore the star-formation activity likely dominates 
the submillimeter emission. It is possible that SSA22-AzTECl is the first example of a protoquasar 
growing at the bottom of the gravitational potential underlying the SSA22 protocluster. 
Subject headings: galaxies: formation — galaxies: starburst — submillimeter: galaxies — infrared: 
galaxies — X-rays: galaxies — quasars: general 



1. INTRODUCTION 

Mounting evidence suggests that submillimeter- 
selected galaxies fSMGs. lSmail et aT1ll997l : lHughes et al.l 
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are the most massive, gas-rich sys- 
2-3 with far-infrared (FIR) lumi- 
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minosity is likely produced by intense episodes of star- 
formation, and therefore it is likely that SMGs are under- 
going rapid growths of their stellar components. A main 
catalyst for starburst i s believed to be m ajor mergers of 
gas-rich galaxies (e.g.. ISmail et al1 ll998). The merger 
event is probably an efficient mechanism to transfer cold 
gas into the nuclear region and fuel the star-formation 
and nuclear acti vity, as suggested by hydrodynamical 
simul ations (e.g., iHopkins et alll2006t iNaravanan et all 
l20lQh . 

An outstanding issue surrounding the SMG popu- 
lation is that dusty starburst galaxies at high red- 
shift can relate to early growth of super-massive black 
holes seen at the present day. It is believed that 
radiation pressure from accretion disks around super- 
massive black holes efficiently strips away the ob- 
scuring gas and dust that fuel star-formation activi- 
ties in their host galaxies, making accreting massive 
black holes an import ant regulator of galaxy forma- 
tion and evolution (e.g ., iSilk fc Reesl 11998 iFabiad 119991 : 
iDi Matteo et al.l 120051) . Early (sub)millimeter surveys 
have revealed that host galaxies of powerful active galac- 
tic nuclei ( AGNs) at high-z are often submillimeter- 
bright (e.g., iMcMahon et all 119941 : iDunlop et all 119941 : 
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stantial fraction (20%-50%) of SMGs have AGNs at 
their center, suggesting a causal connection between 
the SMG phenomena (i. e., bulge formation) and mas- 
sive black hole g r owth (lAlexander et all I2003L l2005ri 
Borvs et all 12004 ISwinbank et al l 12004 iTakata et al 
2006; iMenendez-Delmestre et all l2009TlHainline et al 



2009). SMGs have a volume number density and a 
redshift distribution similar to those of quasars (e.g., 
IChapman et al.l 120051 ). allowing many to propose a hy- 
pothesis that the populations of SMGs and quasars can 
be evolutionarily linked. Moreover, black hole masses 
for SMGs appear to be systematically sm aller than those 
found in quasars ([Alexander et al.l 120081 ). implying that 
some fraction of SMGs with obscured X-ray sources are 
at a protoquasar phase, where the black hole is grow- 
ing more rapidly than in typical galaxies and is about to 
blow off the thick surrounding gas clouds by accretion- 
related outflows from the bl a ck ho l e (e.g . . I Sanders et al l 
1981 IKawakatu et all l200l |200l [20071 iGranato et al l 
2004 120061: lAlexander et all I2005at IBorvs et all 120051: 
Kawakatu fe Wadall2009l )~ 



Theoretical studies of cosmic struct ure formation in a 
cold dark matter (CDM) universe (e.g.. lKauffma^rr][T996l) 
predict accelerated and correlated growths of galaxies 
and black holes in high density environments, such as 
protoclusters of galaxies, which are embedded within the 
most massive dark matter halos collapsing in the early 
Universe. It is natural to expect that the merger rate 
should be enhanced in high-density regions of galaxies, 
and in fact, number excesses of SMGs have been tenta- 
tively claimed towards known overdensities around high- 
redsh ift radio galaxies (jlvison et all 120001 : IStevens et all 
2003). Hence the idea naturally arises that SMGs that 
host accreting massive black holes tightly relate to high-z 
large-scale structures. 

The z = 3.09 SSA22 protocluster was originally iden- 
tified as a significant co ncentration, in reds hift space, 
of Lyman break galaxies (iStei deTet al.l[l9 98). The sur- 
face density of those Lyman break galaxies is ~6 times 
higher than in the field, making this concentration one 
of the highest density regions know to date. Theoreti- 
cal modeling indicates that the protocluster will evolve 
into a rich clust er with a total ma ss of > 10 15 M Q at 
the present day (jSteidel et all 1199 8). The protocluster 
has also been found to contain a factor of ~3-6 overden- 
sityof Lya emitters (jSteidel et al.ll2000l : lHavashino et all 
12004 ). Subsequent spectroscopic follow-up of these Lya 
emitters, which are believed to be star-forming build- 
ing blocks of galaxies, has revealed a large-scale fila- 
menta ry structure extendi ng across ~60 h~ x Mpc or 
more (Matsud a et al.H20"05l ). In this struct ure, many ex- 
tended Lyq-emitting object s (Lyq blobs. iSteidel et all 
[20001: lHavashino eFal|[200l iMatsuda et afl T2004) have 
been found, which are tho ught to be a formatio n 
site of massive galaxi es (e.g.. IGeach et al.ll200l [2007h . 
lUchimoto et al.l (|2008j ) have suggested a density enhance- 
ment of an other massive population of distant red galax - 
ies (DRGs. iFranx et all 120031: Ivan Dokkum et al.ll2003l ). 
Furthermore. iLehmer et all (|2009al) have found a re- 
markable enhancement in AGN in SSA22 compared to 
z w 3 field galaxies. Thus, the SSA22 protocluster is an 



ideal site for studying the co-evolution of massive galax- 
ies and powerful AGNs in a high-density environment. 

Recently a larg e- area 1.1-mm survey toward SSA22 has 
been performed dTamura et al| [2009) using the AzTEC 
1.1-mm camera (IWilson et all I2008al) mo unted on the 
ASTE 10-m telescope ([Ezawa et all 12004 . They have 
mapped a «390 arcmin 2 region of the protocluster and 
revealed 30 SMGs. Among them, the brightest 1.1- 
mm source SSA22-AzTECl (S Ilmm = 8.4t°;g mJy) was 
found near the core of the protocluster, which hints that 
it could evolve into a massive elliptical in this descendant 
of the protocluster. 

However, the multiwavelength counterparts to SSA22- 
AzTECl cannot be easily identified because of the coarse 
beam (30" in a full width at half maximum, FWHM) 
of the AzTEC/ASTE instrument and multiple candi- 
dates within the error circle. For SMG astrometry, cen- 
timeter radio interferometers have often been employed 
so far, but they are not sensitive to SMGs at z > 3 
since radio flux rapidly dims with increasing redshift. 
In contrast, imaging of the dust emission in submil- 
limeter waveband benefits from the strong negative K- 
correction, and its flux density is almost constant for a 
;alaxy with a fixed FIR luminosity at redshifts z « 1-10 
Blain k. Longairl |1993 | ). Recent high spatial resolution 
observations (e.g.JIono et alll2006al : lYounger et al.1120071 . 
120091: iWang et all 120071 I2009D using the Submillimeter 
Array (SMA. IHo et al.ll2004 have uncovered a substan- 
tial population of high-z SMGs that are undetected in 
the radio and optical, suggesting that we would miss 
some high-z SMGs by the radio selection method. An- 
other potential hazard of using wavelengths other than 
(sub)millimeter is that a simple extraction of the nearest 
radio or mid-infrared (MIR) source within an error circle 
could lead to misidenti fication of counterparts to SMGs 
(e.g,. ICowie et al.ll2009D though it is difficult to estimate 
an accurate percentage of misidcntifications. 

In this paper, we present the results from SMA 860- 
/jm observations and multiwavelength properties of the 
brightest 1.1-mm source SS A22- AzTEC 1 found in the 
SSA22 protocluster field. We use spectral energy dis- 
tributions (SEDs) to constrain its photometric redshift. 
We also model an SED in the near-infrared (NIR) to MIR 
to estimate the physical properties of the stellar compo- 
nent. We then utilize X-ray data to identify AGN ac- 
tivity. The multiwavelength study of the 1.1-mm source 
provides a fairly unique opportunity to investigate the 
growth of cosmic hierarchical structures over quite differ- 
ent spatial scales from AGNs (^1 pc) to large-scale struc- 
tures (^10 Mpc): co-evolution of super-massive black 
hole and its host galaxy, and preferential growth of a 
massive galaxy with a powerful AGN in a high-density 
environment in the early Universe, which is one of the 
main concerns in recent astrophysics. 

Throughout this paper, we will assume an £2m = 0.3, 
SIa = 0.7 cosmology with Ha = 70 km s" 1 Mpc" 1 as 
cosmological parameters. An angular scale of 1" corre- 
sponds to 7.6 kpc physical size at z = 3.1. 

2. DATA AND OBSERVATIONS 

2.1. SMA Observations and Data Reduction 

The SMA is an interferometer operating at submil- 
limeter wavebands (230, 350, and 690 GHz) that con- 
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sists of eight 6-meter antennas located at Mauna Kea, 
Hawaii. We used the 'Compact (C)' configuration of 
the SMA. Seven of the 8 antennas were operational in 
the observing track, which gave a range of projected 
baseline lengths of 6-60 meters, which corresponds to 
7-90 kA (i.e., projected antenna separation lengths mea- 
sured in units of an observing wavelength A, which de- 
termine the resultant spatial resolution of interfcromet- 
ric images). The heterodyne receivers equipped with 
a superconductor-insulator-superconductor (SIS) mixer 
were tuned to 355.843 GHz (upper sideband, hereafter 
USB) and 345.843 GHz (lower sideband, hereafter LSB) 
for observing the continuum emission. We obtained 
2 GHz total bandwidth in each sideband. The FWHM 
of the primary bearrF^I at the frequency is 34", which is 
comparable to the AzTEC/ASTE spatial resolution and 
enough to cover the 2a error circle of the AzTEC / ASTE 
positional uncertainty (< 20"). 

The data were obtained on September 6, 2008. The 
conditions were 'excellent' (zenith opacities at 225 GHz, 
T"225GHz = 0.04-0.06, and the typical root-mean-square 
(r.m.s.) phase fluctuation at 350 GHz, A<fi ~ 10°), and 
the double sideband system noise temperatures through- 
out the track were 150-300 K. The phase tracking center 
for the target source SSA22-AzTECl was set to R.A. 
(J2000) = 22 h 17 m 32M, Deck (J2000) = +0°17'35'.'47, 
which is the centroid derived from AzTEC observations 
of SSA22-AzTECl. Each visibility was integrated for 
30 sec, which is short enough to avoid decorrelation 
caused by atmospheric phase fluctuations. The total 
integration time on SSA22-AzTECl was 8.8 ksec or 
2.5 hours. The target source benefits from having a very 
good nearby gain calibrator, a 2-Jy radio-loud quasar 
3C446 (5.6° away from the target). 3C446 was observed 
every 20 minutes to calibrate the gain variation. 3C454.3 
and Uranus were observed at the end of the track for 
bandpass and flux calibration, respectively. The total 
integration time on 3C454.3 was 60 min. 

A baseline error of ~ 0.1 A will yield a systematic po- 
sitional error of an order 0'.'02, which is much below the 
statistical error (0'/13; see $3}. To further check our 
astrometry, a fainter, nearer radio-loud quasar J2218- 
035 (R.A. (J2000) = 22 h 18 m 52?0377, Deck (J2000) = 
-3°35'36"879, which is 3.8° away from SSA22) with a 
precise known position was observed for verification of 
the accurate transfer of the phase solutions and checking 
the astrometric precision. 

The visibility data were calibrated using the iDL-based 
SMA standard reduction package, MIR. We did not need 
to flag out any of the data because of the excellent ob- 
serving conditions. The flux density of a gain calibrator 
3C446 at 860 fj,m was estimated to be 2.12 ± 0.11 Jy. 
The calibrated visibility data in USB and LSB were 
compiled and imaged (Fou rier-transformed) using the 
Miriad (jSault et akl I1995D task, invert. The resul- 
tant size of the natural-weighted synthesized beam is 
3"43 x T/92 (position angle P. A. = 34?3). Since the spa- 
tial frequency coverage is fairly poor because of short ob- 
serving time and the declination very close to zero, the 
resultant map suffers from high side- lobe level of 60%. 

16 The primary beam is a point response function of a single 
dish element of interferometers, and provides a field of view of an 
interferometric synthesized image. 



We extensively CLEAN-ed the dirty image blindly (i.e., 
without constraining the region of interest to model the 
intensity field), restored the clean model, and convolved 
it with the clean beam (a 2-dimensional Gaussian with 
3"43 x l'/92, P. A. = 34?3). The uncertainty of absolute 
flux scaling is estimated to be better than 15%. 

2.2. Multiwavelength Data 

We have multiwavelength data from the radio to the 
X-ray, which were already taken and cover the position 
of SSA22-AzTECl. Here we provide a brief summary of 
the existing data and, in part, newly taken data. 

We use 1.4 GHz (20 cm) radio data o btained towards 
SSA22 that were already published in IChapman et akl 
(200$. The SSA22 field was observed in the B configura- 
tion of the VLA, yielding a ~ 5" synthesized beam. The 
on-source time iinteg = 12 hr resulted in an r.m.s. noise 
level over the AzTEC coverage to be ~ 12 /iJy beam -1 . 

We retrieve MIR photometry data from the Spitzer 
archive to search for a counterpart at MIR wavelengths 
and measure the flux densities. We use 3.6, 4.5, 5.8, 
8.0 /im ima ges obtained with the Infrared Array Cam- 
era fIRAC. iFazio et akl 12001 and 24, 160 [im images 
obtained with t he Multi-band Im aging Photometer for 
Spitzer (MIPS. iRieke et al.ll200l . There is no 70 
data at the position of SSA22-AzTECl. The basic cal- 
ibrated data (BCD) of IRAC and MIPS are processed 
through masking, flat fielding, background matching, 
and mosaicing using the mopex (Mosaicking and Point 
Source Extraction) software, which is a package devel- 
oped at the Spitzer Science Center for astronomical im- 
age processing, along with calibrated data retrieved from 
the Spitzer archive. Then point sources at 3.6, 4.5, 5.8, 
8.0, and 24 /xm are extracted using a mopex pipeline, 
apex, by fitting the point response function (PRF) for 
each band to each source candidate. 

The deep NIR imaging obse rvations of the SSA22 
region have been performed by lUchimoto et akl ( 2008[ ) 
using M0IRCS3 (|Ichikawa et akl [20061 ISuzuki et akl 
2008) mounted on the Subaru telescope, Mauna Kea, 
Hawaii. A part of the NIR data are published elsewhere 
(jUchimoto et al.ll2008h . but we append NIR data that 
are newly taken with MOIRCS/Subaru (Uchimoto et al. 
2010, in preparation). Here we describe a brief summary 
of the MOIRCS observations. The total exposure times 
of the J, H, Xs-band imaging towards SS A22- AzTEC 1 
were 5820, 2765, and 4541 sec, respectively. The stel- 
lar image sizes of the JHKs images are 0"4-0'/5, and 
the limiting magnitudes of the images are J = 25.42, 
H = 24.99, and K s = 25.08 mag (2a) in AB system in 
a 1" diameter aperture. 

In the photometry, IRAC and MOIRCS flux densities 
were measured as follows. In order to achieve higher 
signal-to-noise ratio, we first smoothed the IRAC and 
MOIRCS images so that their PRFs match that of 8.0- 
/jm image. Then we measured the flux densities of the 
NIR-MIR counterparts using a relatively small (3'.'0 in 
diameter) aperture to avoid source confusion. The aper- 
ture corrections were performed for all of the NIR-MIR 

17 Multi-Object InfraRed Camera and Spectrograph devel- 
oped for the 8.2 m Subaru telescope. See also a w eb page at 
http://www.naoj.org/Observing/Instruments/MOIRCS/ for more 
details. 
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data by multiplying a scaling factor. We estimated this 
scaling factor by comparing the 3"-aperture flux and a 
total flux, both of which were obtained for the 8.0-/zm 
image since the 8.0-/im image was less susceptible to 
confusion noise than the other (smoothed) images (see 
Figure |3]). The 8.0-/i, m total flux was measu red with 
SExtractor MAG_AUT0 (jBertin k, ArnoutslH996l ). 

Deep optical images in B, V, R, i\ z' and NB497 
bands taken with Suprime-Cam dMivazaki et al.l I2002T) 
on Subaru is available at the position of SSA22-AzTECl. 
NB497 is a narrow band filter centered at 497.7 nm, 
which is designed to search for strong Lya emission 
in galaxies at z — 3.06-3. 12. The details of th e ob- 
servations are described in Havashi no et al.l (|2004j ) and 
IMatsuda et all (|200l 12005)). The' average profile of a 
point source in the final images has an FWHM of l'.'O. 
All of the optical photometry are measured with a 2'.'0 di- 
ameter aperture, and are corrected for the Galactic red- 
dening of E{B — V) = 0.062. The limiting magnitudes 
after the correction are 28.2, 27.9, 28.3, 28.4, 28.1, and 
27.4 (AB, la) in B, NB497, V, R, i', and ^'-band, re- 
spectively. 

Rec ently, iLehmer et all (|2009allbl ) and iGeach et afl 
(2009) have reported the ultra-deep (400 ksec) Chan- 
dra/ ACIST observations towards SSA22 (Chandra Deep 
Protocluster Survey) . The survey covers a solid angle of 
« 330 arcmin 2 , and revealed 297 X-ray sources, which 
are listed in the "Main Source Catalog". The typical 
sensitivity limits were approximately 4.8 x 10~ and 
2.7 x 10 -16 erg cm -2 s -1 for the 0.5-2 keV and 2-8 keV 
bands, respectively. These limits correspond to rest- 
frame 2-8 keV and 8-32 keV luminosities of 3.7 x 10 42 
and 2.1 x 10 43 ergs s _1 , respectively. 

3. RESULTS AND ANALYSES 

3.1. SMA Results 

In Figure [TJ we found a clear (8.4cr) emission fea- 
ture (AzTEC J221732.42+001744.0; we hereafter call 
the SMA source SSA22-AzTECl for convenience) lo- 
cated ~ 8" northward of the phase tracking center, 
strongly suggesting that this is the true counterpart 
to SSA22-AzTECl. The resultant r.m.s. noise level is 
1.4 mjy beam -1 , which is derived from a region away 
from the 860 /im counterpart. The peak intensity at 
860 /im is 11.8 mjy beam -1 . The position is R.A. 
(J2000) = 22 h 17 m 32?42, Deck (J2000) = +0 o 17'44'/01. 
The inset panel of Figure Q] (lower-right) shows the 
860 /im continuum image obtained towards the test 
source J2218-035, a strong radio quasar with an accu- 
rately known position, which was clearly detected with 
~ lOcr. The offset of the nominal peak position from 
the phase tracking center is (AR.A., ADecl.) = (O'.'OO, 
— 0'.'31). The statistical positional error is estimated to be 

0'.'13 following the formula A0 sta t - v^maj#min/2 SNR, 
where # ma j and 8 min are the FWHM along the major 
and minor axis of the synthesized beam, respectively, 
and SNR is the signal-to-noise ratio of the source in the 
cleaned map. These suggest that the phase transfer is 
applied correctly and the absolute positional uncertainty 
through this SMA observation is estimated to be < 0'.'33. 

We then investigate the flux density and the spatial 
extent of the source more accurately in the wy-domain. 
Figure [2] shows the visibility amplitudes versus projected 



baseline length for SSA22-AzTECl, along with those ex- 
pected for axisymmetric 2-dimensional Gaussians with 
FWHM = 0'.'5, l'.'O, and 2'.'0. The amplitudes are well 
described by a constant value (S l 860/xm = 12.2 ± 2.7 mjy) 
as a function of projected baseline, suggesting that the 
source is not spatially resolved with the SMA 2" beam 
and the upper limit on the source size is < 1". The 
flux densities estimated from the map and the visibilities 
are slightly different (and still consistent) even though 
the source is not spatially resolved. The main difference 
is that the image is synthesized from the vector sum of 
all measurements of the visibilities while the total flux 
estimate in the spatial frequency domain is a scalar esti- 
mate of the vector-averaged subsets of the visibility data, 
which is prone to noise-boosting by the fact that the 
visibility amplitude is always a pos itive quantity (e.g., 
iThompson. Moran. fc Swensonll200ll ). 

Considering the flux density at 1.1 mm is <Si.i mm = 
8-415'q m Jy> which is corrected for flux boosting, the 
860-to-1100 /im flux density ratio is estimated to be 
5 , 860/im/S'i.i mm = 1.45lg45. We scale our result to 
compare the flux ratio for 890 /im SMA follow-up of 
AzTEC sources, and find S$go^ m / Si.i mnl = 1.4 ± 0.4 
for SSA22-AzTECl. This value is consistent with 1.6 ± 
0.7 for the SMA-identified AzTEC sources reported by 
lYounger et all ([2007I . [2009h 

3.2. Possibility of Gravitational Lensing 

The intrinsic flux density of SSA22-AzTECl may be 
lower if there is an amplification by gravi tational lens- 
ing. While recent theoretical works (e.g.. iPaciga et all 
2009) imply that the fraction of «10-mJy submm 
sources magnified by a factor of > 2 is relatively small 
(~5%), there are many SMGs that are e ither known 
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small spatial size of SSA22-AzTECl measured with SMA 
(Figure [5]) suggests that the high brightness of SSA22- 
AzTECl is unlikely due to a strong gravitational lensing 
effect caused by a foreground massive galaxy that would 
happen to be aligned in the direct line of sight. More- 
over, if not in the direct line of sight, a lensing object 
must be found at an angular separation comparable to 
the Einstein radius. If a source at z > 1 is located close 
to a foreground lensing object at z ~ 0.3, which is the 
favored geometry for efficiently lensing the background 
source, Einstein radii are ~ l"-2" and ~ 30"-40" when 
the lensing object is a galaxy (dark halo mass of 10 12 M Q ) 
and a cluster (10 15 M Q ), respectively. But, we do not see 
such a low-z massive galaxy ~1"— 2" away from SSA22- 
AzTECl although it would be easily detected in the Sub- 
aru Ks image. In SSA22 we find neither X-ray clusters, 
which are often found at z ~ 0.3, nor higher- z [z ~ 1) 
Sunyaev-Zel'dovich clusters, which should be observed in 
the AzTEC 1.1 mm map if exists. It is also unlikely that 
the SSA22 protocluster at z = 3.1 itself is a lensing ob- 
ject because the SSA22 protocluster is by definition far 
from being virialized and should not have deep poten- 
tial well similar to those of low-z clusters. It is therefore 
likely that SSA22-AzTECl is not strongly magnified by 
gravitational lensing. 
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3.3. Multiwavelength Counterparts 

At the position of the 860-/zm source, we find a sig- 
nificant counterpart in the radio to MIR as well as in 
the X-ray, but the 860-/zm source has no counterpart in 
the deep Subaru optical images. The flux densities and 
upper limits at radio to optical wavelengths are given in 
Table [TJ The postage stamp images at radio (20 cm) to 
near-ultraviolet (230 nm) wavelengths are presented in 
Figure OH 

Although there are multiple counterpart candidates 
in the radio map, we identify a radio counterpart 
(<Si.4GHz = 42 ± 12 /zJy) at the position of the 860-^m 
source. The spectral index between 860 (350 GHz) 
and 20 cm (1.4 GHz) is a? 5 4 ° = 1.03 ± 0.15, suggesting a 
redshift of z ~ 3-4 on the basis of the r edshift-estimator 
proposed bv lCarilli fc Yunl (fl999l l2000h . 

We find counterparts to the SMA position in all IRAC 
bands and MIPS 24 fim. However, we see no significant 
counterpart at 160 mainly due to heavy blending 
from infrared sources, which are located ~ 15" south- 
ward from the 860- /urn position (see Figure [3]). These 
infrared sources have radio counterparts, and are the 
brightest at 24 ^m in the SMA field of view. SSA22- 
AzTECl has red color in the Spitzer/IRAC bands, which 
is consistent with a z ~ 3 SMG. Figure 2] shows the 
color-color (IRAC 5.8-to-3.6 and 8.0-to-4.5 /im flux ratio) 
diagram for this source, SMGs with millimeter carbon 
monoxide (CO) or MIR Spitzer/IRS spectroscopy, and 
infrared-lumi nous quasars from the Spitzer First Look 
Survey (FLS. lLacv et~al|[200l . The 5.8-to-3.6 and 8.0- 
to-4.5 fim flux ratios of this object are log (S^.s/S^.q) = 



-0.03 



and log (Ss.o/Sa.5 



0.351q'o3, respectively. 



0.46; 

Indeed, it is found near the transition region between 
z < 3 and z > 3 SMGs with secure redshifts. It is seen 
in Figure [4] that infrared-luminous quasars identified in 
the FLS survey also have red IRAC colors, and it seems 
the z > 3 SMGs indeed have colors that overlap with the 
infrared-luminous AGNs with a power-law spectrum in 
MIR (thick magenta line). We note that SSA22-AzTECl 
has the IRAC colors consistent with a z ~ 3 SMG, but 
all z > 3 SMGs also have IRAC colors very similar to 
the infrared power-law AGNs. 

The SED of SSA22-AzTECl is well constrained at 
wavelengths > 3 /im, but drops out at a wavelength 
of 1 /im (J > 24.6 in AB, 2a) or shorter. We ten- 
tatively detect counterparts to SSA22-AzTECl in H- 
(3a) and ATg-bands (4er), whereas we see no emission 
in J-band. The observed color of SSA22-AzTECl in the 
NIR bands (J — K > 1.0) may meet the color criterion 
(J - K > 1.4 in AB) for distant red gala xies (DRG, 
iFranx et all 120031: Ivan Dokkum et all [2003h . which are 
often found in the Universe at z ~ 2 and thought to 
be o ne of the most massiv e populations at this epoch 
(e.g.. [Kaiisaw a et al.l [2006). This extremely red spec- 
trum in the NIR-MIR strongly suggests that SSA22- 
AzTECl can have a large obscuring column in front of 
the stellar component and/or an AGN. 

We find no significant emission in the Subaru opti- 
cal B, V, R, i', z' and NB497 images. The photo- 
metric data from the Subaru observations are tabulated 
in Table [TJ The non-detection in NB497 does not im- 
mediately imply that the SMA counterpart is not at 
z « 3.1 since, due to significant extinction, Lya emis- 



sion in ~ 50% of SMGs with a bright optical coun- 
terpart would not be stron g enough to be detect ed in 
narrow-band imaging (e.g.. IChapman et al.l I2005T ). In 
addition, there are no counterparts in Hubble optical 
(HST/ACS F814W) and GALEX near/far-ultraviolet 
images (see Figure [3]). The properties in the rest- 
frame optical images are very similar to those of submil- 
limeter sources with secure (sub)millimcter i nterferome- 
ter identifications GOODS 850-5 (or GN10. iPope et all 
[2001 IWang et all [20071 [20091: IDannerbauer et al.l | 200£ ' 
iDaddi et all I2009al ). HDFN850TT~ lHughes et all 11995 



Do wnes et all 11999b Dunlop et al.l 



2009), SXDF850.6 (|Hatsukade et al 
identified AzTEC sources in the 



2004; ICowie et all 
12010I) . and SMA- 
COSMOS field 
(|Younger et all [20071 [20091) . 

It is worth noting that in the proximity of SSA22- 
AzTECl (~ 10"), we find a large emission nebula in 
the Subaru/NB497 image, which is a corroborated can- 
didate for a Lya blo b at redshift z = 3.1 (LAB36, 
lHavashino et alll2004|) . In Figure[TJ we show a false-color 
optical (B, NB4 9 7, V) image around SSA22-AzTECl. 
lHavashino et all ((2004) have found 74 Lya blobs over 
their surveyed area in SSA22 (699 arcmin 2 ), 59 of which 
are concentrated in the so-called high-density region 
(302 arcmin 2 ) of SSA22, in which SS A22- AzTEC 1 is lo- 
cated. The probability that an object has a chance to 
be closely (< 10") associated with one of the 59 Lya 
blobs over a 302 arcmin 2 region is 1.7%. The separation 
angle of ~ 10" corresponds to a projected separation of 
~ 80/i^q 1 kpc in proper scale at z = 3.1. 

We find a significant X-ray counterpart (~ 20 counts), 
which is listed as the source #120 in the Main Source 
Catalog of the Chandra Protocluster Survey, coincident 
with the SMA source position, suggesting the existence of 
an accreting massive black hole buried in the dust clouds. 
The position of the X-ray source is R.A. (J2000) = 
22 h 17 m 32?42 and Deck (J2000) = +0°17'43'/9, and its as- 
sociated positional error is 0"54, which is estimated from 
the 80% confidence interval. The X-ray counterpart has a 
hard X-ray spectrum (the band raticF 5 !. BR — 5.5S^ 2 ' 2 q, 
or the effective photon indeaP^I. r o ff = — 0.34+q g£) and is 
quite bright in the 2-8 keV band (but faint in the 0.5-2 
keV band; see also Table [TJ). In fact the X-ray spectrum 
is harder than that found for an AGN in a typical SMG 
among a sample of X-ray lumi nous SMGs at z = 0.6-2.9 
([Alexander et all [2003. 2005b), suggesting that it hosts 
one of the most heavily obscured distant AGNs known. 
We will further investigate the X-ray properties in § 13.61 



3.4. Photometric Redshift Estimates 

The important next step is to estimate the redshift 
of SSA22-AzTECl. It is extremely difficult to estimate 
its redshift through conventional optical or NIR spec- 
troscopy because of the faintness in the optical and NIR. 
Here we estimate the photometric redshift (photo- z) by 
fitting the photometric data points to SED templates 
from stellar population synthesis models (§ 13.4. 1|) and li- 
braries that simulate MIR-to-radio emission from warm 
and cold interstellar medium (§ 13.4.21) . 

18 Ratio of the count rates in the 2.0-8.0 and 0.5-2.0 keV bands. 

19 A parameter for a power-law spectrum defined such that / oc 
v — reft, where / and v are flux and frequency, respectively. 
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3.4.1. Photometric redshift from NIR data 

We use synthetic SED models of the composite stel- 
lar populations (CSP) with e xponentially-declining star - 
formation activities from the iBruzual fc Charlotl (|2003f ) 
library. We use the models to fit the IRAC and MOIRCS 
photometry as well as the 2a upper limits in i' and z' 
bands listed in TableQ] and compute the reduced-x 2 (i.e., 
X 2 /dof, where dof is the degree of freedom) as a function 
of redshift, 



1 



fgjobs gmodcl^2 



N-p 



(1) 



where N is the number of photometric data points, p 
is the number of free parameters, S° h ? and 5'™ odcl are 
the flux densities from observation and templates, re- 
spectively, and (Ti is the uncertainty associated with each 
photometric datum. Note that we do not use the pho- 
tometry data at wavelengths A Q b s > 3(1 + z) /im since 
emission from non-stellar components, such as hot dust, 
may conta minate th e flux densities at A rcst > 3 /im. We 
assume the lSalpeterl (|1955f ) initial mass function and the 
solar metallicity (Z = 12©). We treat redshift, stel- 
lar mass, age, star-formation timescalj 2 ^. and extinction 
as fre e parameters. We make use of the ICalzetti et al.l 
(2000) extinction law to account for the extinction with 
color excess of E(B — V) — 0.0-3.0, which corresponds to 
visual extinction of Ay » 0-12. We consider the redshift 
range from z — to 7. 

The photo- z derived from fitting the stellar SED of 
SSA22-AzTECl is poorly constrained. We found a min- 
imum of rcduccd-x 2 of 0.3 at z p hot — 3.9 with the con- 
fidence interval of z > 1.0 (99% confidence level) or 
z > 2.8 (68% confidence level). This means that the stel- 
lar population synthesis models at any redshift z > 2.8 
can easily reproduce the MOIRCS and IRAC SED of 
SSA22-AzTECl and we virtually can not place strong 
constraints on the photometric redshift from fitting the 
MOIRCS and IRAC data to the stellar models. 

The primary reason for the poor constraint is that 
SSA22-AzTECl has no notable features in the rest-frame 
ultraviolet-to-NIR, such as the Lyman break, Balmer 
break, or the 1.6-/im stellar bump, which play critical 
roles in the ultraviolet-to-NIR photometric redshift es- 
timates. As mentioned in the previous section, SSA22- 
AzTECl has an extremely red color in the rest- frame 
ultraviolet-to-NIR, suggesting a large extinction that 
makes a steeply declining SED towards shorter wave- 
lengths. Such a steeply declining SED could cause errors 
in photometric redshift. 

An example of such an error can be found in a series 
of papers on an SMC, GOODS 850-5 (or GN10). Wang 
et al. (2007, 2009) have reported extremely faint SED of 
this object over the rest-frame ultraviolet-to-NIR bands 
(an extremely deep limit to A"s(AB, 2a) < 27.2 mag or 
46 nJy). They also found that the photo-z estimates for 
the ultraviolet-to-NIR data show z p hot ~ 6.5, whereas 
the photo-z from its FIR-to -radio (24 /im-20 cm ) SED 
shows Zp hot ~ 4. Eventually. iDaddi et al.l (|2009b[ ) blindly 

20 CSPs are modeled assuming exponentially declined star- 
formation activities. The star-formation timescale is defined as 
e-folding time-scale of the star-formation activities. 



searched for CO in this object and found the spectro- 
scopic redshift of zqo = 4.042, which as a result rejected 
the photo-z estimated from the rest-frame ultraviolet-to- 
NIR. 

3.4.2. Photometric redshift from MIR to radio data 

Another way to constrain the redshift is to use the FIR- 
to-radio SED. The ratio between MIR-to-(sub)mm and 
radio flux densities depends on the redshift. This is be- 
cause in the local Universe the radio luminosity is tightly 
correla ted with the FIR l uminosity (FIR-to-radio corre- 
lation iHelou et all 119851 : iCondonl 099% IHelou fc Bicavl 
119931: lYun et al.l 120011 ) a nd this r elation seems to hold 
at h i gh redshifts (e.g.. [Garrettl 120021: lAppleton et al 
2004HIbar et al.l 120081: ISevmour et alj 12009: llvison et al 



2010®" There have been a number of studies exploit- 



ing this relatio n to constrain the photometric redshift 
of SMGs (e.g.. ICarilli fc Yunl Il99l IgOOCt IHuehes et all 
200llAretxaga et al.1120031 120051 120071 IWang et alJI2009f 
Daddi et all 12009 a*). and this method can be efficient 



when a galaxy suffers from extremely heavy extinction 
in the rest-optical/NIR bands. 

We use the available photometry at 24, 860, 1100 /im, 
and 20 cm data, and e mploy the SED library from 
IMicha lowski ct al. (2010]) to estimate the redshift of 
SSA22-AzTECl. The library consists of SED mod- 
els, which were developed using the GRASIL code 
(jSilva et al.l 119981 ). th at were being fit to t he photomet- 
ric data points (e.g., Mainline et al.ll2009fl of 76 SMGs 
with spectroscopic redshifts (|Chapman et al.| [2005) . The 
mean FIR luminosity, stellar mass, and redshift, aver- 
aged over the 76 SMGs, are L F ir = 1O 127 L , M star = 
10 11j Mq, and z = 2.0, respectively. To eliminate 
poorly constrained SEDs from the library, we chose SEDs 
with robust 24- /jm photometry (i.e., SEDs that are well- 
constrained at MIR wavelengths and the Wien tail of 
the FIR dust emission). We do not use the SEDs that 
strongly underpredict flux densities at MIR or submil- 
limeter wavelengths (A bs = 3.6-8.0, 850 /im). For our 
final SED templates, we employ 35 SEDs and an aver- 
aged SED that is computed by averaging the 35 SEDs. 
The mean FIR luminosity, stellar mass, and redshift of 
the 35 SMGs are L Fm = 10 12 S L Q , M star = 10 116 M Q , 
and z = 2.0, respectively, which are consistent with those 
averaged over the 76 spectroscopically-identified SMGs. 
Following Equation[TJ we calculate the x 2 values between 
the photometric data of SSA22-AzTECl and each SED 
model with varying redshift and luminosity as free pa- 
rameters. The x 2 n t to the data at the observed wave- 
lengths of A bs > 24 /jm provides a good fit, but the 
X 2 fit is poor when we include the A b s < 8 /im data. 
The A bs < 8 /im data likely suffer from heavy extinc- 
tion by dust if the data at A b s < 8 /im represent the 
stellar component of SSA22-AzTECl. Furthermore the 
SED models are biased agai nst an optically faint pop - 
ulation of SMGs as noted by IMichalowski et al.l (2010). 
Therefore we use the photometric data at A b s > 24 /im 
hereafter. 

Figure [SJi (top) shows the reduced x 2 values for the 
35 SED models with 24-/im data (thin curves) and a 
model averaged over the 35 SEDs (a thick curve), as a 
function of redshift. Figure [S^, (bottom) shows the best- 
fit FIR luminosity that gives the least x 2 value at each 
redshift. The x 01 the bulk of the SEDs have their 
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local minima at z « 3. The x 2 minima, however, scatter 
broadly in redshift space probably because we have only 
tw o degrees of freedom and many of the SEDs presented 
in Michalowski et al. (2010) remain poorly constrained 
especially in the MIR band. 

AGN-dominant galaxies could have high dust tempera- 
tures, which increases the MIR luminosities and shifts the 
FIR peak of SEDs to shorter wavelengths. Thus redshift 
estimates using such SEDs may tend to be biased toward 
higher-z. In fact, there are just three SED templatef! 
having their local minima at z « 3.8. Independent 
measurements of the ir spectra at rest-frame ultraviolet 
to M I R wavelengths dSwinbank et alj|2004: ITakata et~aL 
2006; Menendez-Delmestre et al. 2007 n2009[1PoDe et al.l 
20081 : Mainline et al.l 120091 ) have suggested that they are 
probably AGN-dominant (i.e., AGN is responsible for 
>50% of the bolometric luminosity of its host). The 
AGN of SSA22-AzTECl does not appear to dominate 
the host bolometric luminosity as discussed in the fol- 
lowing section (§ I3.6|) . suggesting that the fits to the 
AGN-dominated SEDs might not be appropriate. On 
the contrary, SED templates with lower dust tempera- 
ture can decrease the best-fit redshift. In our SED fits, 
local minima are evident at z ~ 1.5. The 9.7-/im sil- 
icate absorption feature falls to the 24- fim. data point 
at z ~ 1.5, slightly improving the x 2_n t f° r this partic- 
ular redshift. However, the x 2 values at local minima 
around low redshifts z ~ 1.5 are much larger than those 
for SEDs with best-fit redshift of z ~ 3-4, and the stellar 
SED fits do not favor redshifts z ~ 1. This suggests that 
the redshift of z ~ 1.5 is unlikely. 

It is thus reasonable to focus on the averaged SED. 
The photometric redshift from the averaged SED is 
Zphot = 3.19^035 (the error bar is estimated from the 
99% confidence interval), consistent with the redshift 
of the protocluster (z = 3.06-3.12). The derived FIR 
luminosity is Lfir = 1-9 jl] g X 10 13 L Q (the error bar 
is estimated from the 99% confidence interval). The 
large FIR luminosity impl ies a star- f ormat ion rate of 
~ 4 x 10 3 M Q yr" 1 based on lKennicuttl (|1998ft . The best- 
fit averaged SED along with the photometric data points 
is shown in Figure [5Jd. We note that if we take the best- 
fit SED model (SMM J123711.98+621325.7) that has the 
least minimum-x 2 among all of the SEDs, the inferred 
photometric redshift is z ~ 3.18^Q 4^, again consistent 
with the protocluster redshift. We will hence assume 
z = 3.1 for the redshift of SSA22-AzTECl in the rest of 
the analyses. 

3.5. NIR-to-MIR Spectral Energy Distribution 

While the iMichalowski et al.l ()2010D SED model is in 
good agreement with the data at observed wavelengths 
of A bs ^ 24 /im, we see a significant dimming/break in 
the SED at A Q b s < 8 [im compared with the SED model. 
This is probably because of (i) extremely heavy dust ex- 
tinction of the stellar component, (ii) low-mass stellar 
component, and/or (iii) a steep power-law spectrum in 
rest-frame NIR bands due to hot dust heated by a buried 
AGN. 

To further investigate the possibilities (i) and (ii) 

21 SMM J105238.30+572435.8, SMM J123600.15+621047.2, 
SMM J163650.43+405734.5 



shown above, we use stellar population SED models to fit 
the MOIRCS and IRAC photometric d ata. T he method 
used here is the same as introduced in § 13.4.11 but we use 
pop ulation synthesis m o dels f rom the IMarastonl ()2005[ ) 
and iBruzual &: Chariot! (j2003| ) libraries and fix the red- 
shift to z = 3.1. 

The best-fit SEP mod els from IMarastonl (|2005f ) and 
IBruzual fc Charlotl (|2003l ) are essentially consistent, and 
each of the best-fit SED suggests Ay = 3.4±0.2 mag, age 
of 2.5± 3 /° x 107 y r > and stellar mass of M star = 7.3i 7 _ 7 x 
10 10 M Q (x 2 — 0.963, from the Maraston model) and 
A v = 3At° Q A 3 mag, age of 2.5±^ x 10 7 yr, and M star = 
8.2±l 5 x 10 10 M Q (x 2 = 0.886, from the Bruzual & Char- 
lot model). The uncertainties represent the la confidence 
interval. We model CSPs with e-folding time-scale of 
star- formation activities ranging tsf = 0.1-20 Gyr, but 
we cannot place a good restriction on the time-scale tsf 
because the best-fit SED model is fairly young compared 
to Tgp. We note that the ages from the stellar model 
fits are known to be somewhat unreliable, especially for 
stella r populations in the early Universe ( Si mple v et al.l 
120051 ). The extinction is comparable to but slightly larger 
than those foun d in SMGs with rest-fr ame optical spec- 
tra obtained bvlSwinbank et alJ (120041. Ay = 3.0 ± 1.0) , 
iBorvs et all (I2005L Ay = 1.7 ± 0.2), ITakata et alJ pOOa 
A v = 2.9 ±0.5), and [Michalowski et al.l (120101 median 
A v = 2.03 ± 0.95). IBorvs et alJ (|2005f ) have studied 
stellar components of 13 SMGs found in the GOODS- 
N field, and found the stellar masses of ~ 1O 11 -1O 12 M 
and a mean stellar mass of ~ 2 x 10 11 M Q . The stellar 
mass of SSA22-AzTECl is slightly smaller than those for 
a typical SMG in GOODS-N. But it is larger than that of 
coe val ultraviolet- s electe d star-forming galaxies reported 
by iShaplev et all (I2TM logM a t ar /M ^ = 10.32 ± 0.51 
for z ~ 2) and iMagdis et all (|2010L log M star /M Q = 
10.62±0.11 for z ~ 3), suggesting that SSA22-AzTECl is 
more massive than a coeval typical star-forming galaxy. 
The SED fit improves (x 2 = 0.106) if we employ a super- 
solar metallicity Z = 2.5Z Q although the best-fit param- 
eters do not change drastically (Ay = 2.4+qJ, age of 

1.6j^6 x 10 7 yr, and M star = 3.3±g 5 x 10 10 M Q ). Such 
higher metallicities are often favored in mod eling SEDs of 
local ultra-luminous infrared galaxies (e.g.. iFarrah et all 
2005). We note, however, that since the number of the 
photometric data points at < 8 ^m is limited to just six, 
the parameters (i.e., Ay, age, and M sta r) that determine 
the stellar SED are highly degenerate and could have 
large systematic uncertainties. 

3.6. Luminous AGN Buried Deeply in SSA22-AzTECl 

Here we investigate the detailed properties of the X- 
ray counterpart referred in § 13.31 We utilize the pimms 
packag<£3to model the X-ray spectrum that accounts for 
the small effective photon index of r c ff = — 0.34^o'gi- 
We simulate the intrinsic power-law spectra with pho- 
ton indices Y = 1.4-2.0 behind hydrogen columns in the 
range A^h = 1 x 10 23 -4 x 10 24 cm" 2 . The spectra are red- 
shifted to z = 3.1 and adjusted for absorption through 
the Galactic hydrogen column towards SSA22-AzTECl 

22 PIMMS (Portable, Interactive Multi-Mission Simula- 
tor) is maintained by Koji Mukai, and is available at 
http : //heasarc .nasa.gov/docs/sof tware/tools/ 
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(4.1 x 10 20 cm- 2 . iKalberla et all 120051) . Other compo- 
nents such as reflection and scattering are not consid- 
ered. The band ratios between the soft (0.5-2 keV) and 
hard band (2-8 keV), BR, are estimated using pimms 
(v. 3.7a) appropriate for the Chandra / ACIS-I Cycle 8 
observations. For BR > 3, band ratios do not sub- 
stantially depend on intrinsic photon index, but depend 
largely on hydrogen column density. If the AGN has 
intrinsically a power-law spectrum with a photon in- 
dex of r = 1.8, which is a typical value for nearby 
AGNs (jTozzi et al.ll2006f l. the hard spectrum of SSA22- 
AzTECl strongly suggests that the hydrogen column has 
JVh ~ 1 x 10 cm -2 and is almost Compton-thick. It 
should be noted that the band ratio works well for AGNs 
with Nn < 5 x 10 23 cm -2 in determining the hydrogen 
column density, but it becomes less reliable for higher 
column densities because other components such as re- 
flection and scattering make the X-ray spectrum more 
complex than just an absorbed power law. However, the 
X-ray spectral slo pe is as flat as t h e heav ily-obscured 
SMGs reported by lAlexander et al.l (|2005bf ) which also 
appear to have Ah ~ 10 24 cm" 2 , based on more com- 
plex X-ray spectral fitting. 

Assuming that the X-ray source is associated with 
SSA22-AzTECl and it is located at z = 3.1, the 
absorption-corrected luminosity in rest- frame 0.5-8 keV, 
Lx, is estimated following the expression, 

i x -4^ 2 ^ bs (l + z) r " 2 , (2) 

where cZl = 26.4 Gpc is the luminosity distance at z = 3.1 
and F-£ bs is the absorption-corrected flux in the observed- 
frame 0.5-8 keV. We find L x w 3 x 10 44 ergs a~ l , 
which is slightly larger than those found in z ~ 2 AGN- 
cl assified SMGs (L x ~ 10 44 er S s s_1 j see Figure 8 
in lAlexander et al.l l2005bl ) . It is therefore likely that 
SSA22-AzTECl hosts a heavily obscured AGN, which 
is one of the most luminous objects found among the 
SMG population. 

Given the large X-ray luminosity, the bolometric out- 
put from the central AGN in SSA22-AzTECl can heat up 
the interstellar dust and it may contribute to a substan- 
tial fraction of the FIR luminosity of the host galaxy. It 
is becoming clear that AGNs are likely present in a signif- 
icant fraction of SMGs but their contrib ution to the FIR 
emissi o n is in many cases minor (e.g., [ Alexande r et all 
20054 iMenendez-Delmestre et al.l 120091 : Mainline et al.l 
20091 iMichalowski et al.l 12010ft . Nevertheless, it is im- 
portant to investigate the origin of the FIR luminosity 
of SSA22-AzTECl (and hence the intrinsic bolometric 
output from the AGN) because the AGN-to-host lumi- 
nosity ratio can be used as an indicator of AGN activity 
and may provide a hint to the evolutionary phase of the 
central massive black hole. We now explore what powers 
the bolometric luminosity of SSA22-AzTECl — massive 
star-formation or an AGN. 

The unabsorbed X-ray to FIR luminosity ratio of 
SSA22-AzTECl is Lx/L Fm = 0-004. This value i s very 
similar to that found by lAlexander et al.l (f2005b) (me- 
dian Lx/I/fir ~ 0.004), but it is an order of magnitude 
smaller than tho s e found in quasars (Lx /L F m, ~ 0.05, 
lElvis et all 119941: lAlexander et all 12005b). The reason 
why the luminosity ratio in SSA22-AzTECl is so much 
smaller than those in quasars is the extreme FIR lumi- 



nosity of SSA22-AzTECl. In order to estimate the bolo- 
metric luminosity of the AGN, we introduce the bolo- 
metric correction, which is defined as /boi = L^™ /Lx, 
where is the bolometric luminosity of the AGN. 

The bolometric correc tion for the X-ray band is known 
to be ~10-100 (e.g., lElvis et al.l 119941 ) and to have a 
mo derate dependence on the bolometric AGN luminos- 
ity. iMarconi et al.l (|2004l ) have investigated AGNs in the 
local Universe, and found an empirical relation between 
L£j N and / bo i in 2-10 keV band as 

log /boi = 1.54 + 0.24/3 + 0.012£ 2 - 0.0015£ 3 , (3) 

where £ is a bolometric luminosity measured in units 
of 1O 12 L and is defined as £ = logLj^f 1 - 12 = 
log /boi + \ogL^ whcV - 12. We correct the SSA22- 
AzTECl's 0.5-8 keV luminosity to 2-10 keV, and ob- 
tain L x ~ 10kcV = IO W - 7 L assuming T = 1.8. We solve 
Equation ([3]) in terms of £ by substituting log /boi = 
£ - log L x _10keV + 12 = £ + 1.3 and find £ = 0.33. 
Wc then have L^ N 10 12 ' 3 L Q and / bo i ~ 40 for 
the AGN of SSA22-AzTECl. The bolometric correc- 
tion is consistent with those measured for high-z AGNs 

(/bol 10kcV » 35. lElvis et al.l H99l. 

This unabsorbed bolometric luminosity of the AGN is 
extremely large, and as luminous as the m ost luminous 
objec ts in the local Universe (e.g., Mrk 231. iSoifer et al.l 
1987). Nonetheless, the bolometric AGN luminosity can 
only account for ~10% of the huge bolometric output 
of SSA22-AzTECl (L Fm « 2 x 10 13 £ s ). If we adopt 
the highest estimate for a bolometric correction of the 
X-ray luminous SMGs (/b~! 10kcV « 70), the percentage 
will increase by a factor of ^1.6. This suggests that the 
large FIR luminosity of SSA22-AzTECl is not mainly 
powered by the AGN, but should be attributed to other 
power sources, likely massive star-formation activities. 

4. DISCUSSIONS 

We have shown that SSA22-AzTECl is a hyper- 
luminous starburst galaxy at high redshift, likely z = 3.1, 
and harbors a luminous X-ray source. Here we discuss 
a possible scenario for the growths of the stellar compo- 
nent and massive black hole of SSA22-AzTECl in the 
protocluster environment. 

4.1. A Progenitor of a Massive Galaxy in a Protocluster 

The inferred redshift and the large stellar mass shown 
above suggest that we are possibly witnessing the vio- 
lent forming stage of a massive early-type galaxy that 
nucleates at the bottom of the gravitational potential 
underlying this protocluster, which can be the archetypi- 
cal formation site of massive galaxies predicted from the 
current standard model of structure formation in a CDM 
universe. There is a claim that the bright-end of the 
evolved systems selected at NIR wavelengths is well pop- 
ulated by z ~ 2 but much less in z ~ 3 protoclusters 
around powerful radio galaxies, implying that the bright- 
end of the evolved, likely massive (~1 n M^) galax- 
ies h ad first appeared between z — 2-3 (Kodama et al. 
120071) . The time duration from z = 3 to 2 is wl.l Gyr 
whereas just 0.01 Gyr (i.e., comparable to a lifetime of 
a massive star) is required for SSA22-AzTECl to have 
10 1:l Mq of stars if the stars could be produced at cur- 
rent star-formation rate (~ 4 x 10 3 M Q yr _1 ) and add 
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another ~(2-3) x 10 10 Af Q to this system. Although we 
have no information on the mass of the gas reservoir in 
SSA22-AzTECl, it is possible for this object to evolve 
into a M st ar ~ lO n M0 system by z = 2. A search for a 
molecular reservoir in this object will allow us to assess 
this possibility. 

4.2. Growth of a Massive Black Hole in a Protocluster 

Comparing the black hole and stellar mass is a good 
method for understanding the evolutionary status of the 
AGN in SSA22-AzTECl. While estimating masses of 
black holes (Mbh) is often fairly difficult, there are ways 
to constrain the black hole mass. One of the most plau- 
sible ways using the available data for SSA22-AzTECl 
is to assume the Eddington ratio, r\ (i.e., ratio of the 
bolometric to Eddington luminosity of the AGN) . When 
77 = 1, the AGN is powered by Eddington- limited accre- 
tion onto the black hole, which virtually provides a lower 
li mit on an estima t e of th e black hole mass. 

I Alexander et al.l (|2008l ) have investigated A^bh and 
rj of spectroscopically-identified SMGs in the Chan- 
dra Deep Field-North. They divided these SMGs into 
two samples: one is a sample which consists of SMGs 
that exhibit optical broad emission lines. The other 
consists of SMGs that host X-ray identified obscured 
AGNs with no apparent broad-line region, which is 
more representative of the overall SMG population than 
the broad-line SMG sample. The average redshifts of 
both samples are consistent with the median redshift 
obtained for ra dio-identified SMGs ((z) mc dian = 2.2, 
IChapman et al.l | 2005D. and t heir s tellar masses are well 
constrained by iBorvs et all (|2005t ). For the broad- line 
SMGs, one can directly estimate Mbh and r\ using the 
well-established virial black-hole ma ss estimator (e.g., 
iWandel et aDll99l iKaspi et al.l [2001 . They found the 
mass ratios Mbh/Mqal of the broad-line SMGs and the 
X-ray-obscured SMGs to be ~ 2 x 10~ 3 and w 2.9 x 10~ 4 , 
respectively. The average mass ratio of the broad-line 
SMGs is higher than that of the X-ray-obscured SMGs 
but lower than those of typical unob scured quasars at 
z ~ 2 (M R h/Moa t, ~(5-10) x 10~ 3 . iPeng et all [20061: 
ICoppin et al.l 12008) . This suggests that the broad- line 
SMGs have more massive, rapidly growing black holes 
than the X-ray-obscured SMGs, and they are possibly at 
a transition phase between X-ray-obscured (i.e., typical) 
SMGs and unobscured quasars. 

Although uncertainties are considerably large, the in- 
ferred mass of the black hole in SSA22-AzTECl is 
MaH.Edd 10 7 M Q for Eddington-limited accretion, 

or M BH ~ 2.5 x 10 s (r?/O.2)- 1 M . The black-hole-to- 
galaxy mass ratio of SSA22-AzTECl, which is estimated 
from the black hole mass divided by the stellar mass, 
is Mbh/Mgal = 3.5 x 10- 3 (t7/0.2)- 1 . The mass ra- 
tio is (1.4-3.5) x 10 -3 if the Eddington ratio has a re- 
alistic value between ij — 0.2-0 .5, which are found t o 
be true in the SMG population ( Alexan der et all 2008). 
This mass ratio is ~5-10 times higher than those found 
for the X-ray-obscured SMGs, but it is below the val- 
ues found in the z ~ 2 unobscured quasars. This large 
discrepancy might be attributed to the underestimation 
of the stellar mass, which is difficult to estimate as dis- 
cussed in § 13.51 If we underest imate the stellar ma ss by 
as large as a factor of 2 (e.g., [Shaplev et aLl I2005T) . the 



mass ratio could still be 2.5-5 times larger than those 
of the X-ray-obscured SMGs. A more likely explanation 
is that SSA22-AzTECl has a more massive black hole 
accreting at a higher rate (i.e., higher rf) than the X-ray- 
obscured SMGs, but it is less massive than black holes 
of the quasars. In contrast, the mass ratio is compa- 
rable to those of the broad-line SMGs, suggesting that 
the black hole in SSA22-AzTECl can be in the same 
evolutionary stage as those in the broad-line SMGs. Al- 
though the constraints are very uncertain, it is hence 
possible that the black hole in SSA22-AzTECl is be- 
ing at a tra nsition phase to a quasar, or a vrotoquasar 
phase (e.g. , iSanders et all | 1988c IKawakatu et al.l 120031 
2001 [20071 iGranato et ai] 12001 l2006t lAlexander et al.l 
2005at IKawakatu fc Wa da 2009), where the black hole is 



growing more rapidly than in typical SMGs and is still 
buried deeply in the central gaseous, dusty 'cocoon'. 

What is the importance of finding a luminous ac- 
creting massive black hole, which is one of the most 
obscured AGN found in SSA22 and likely resides at 
the heart of the brightest SMG in the SSA22 proto- 
cluster? Large-area submillimeter surveys have sug- 
gested that SMGs are related to la r ge-scale structures 
at high redshift (IBlain et al.l [20041 ; IScott et all 1200 
Daddi et al.l l2009allbl : ITamura et all 120091 : iViero et al 



2009tlWeifi et alJ2009b|) even though there is a claim that 



the population of SMGs cannot always be used as tracers 
of the most massive dark halos at high-z considering the 
short time d uration (as short as 10 Myr) of their star - 
burst phase (IChapman et al . 2009; Arav ena et aLl feOlO). 
In particular, ITamura et al.l ( 2009 ) have found an an- 
gular correlation between the positions of the bright 
SMGs and the Lya emitters in SSA22, strongly sug- 
gesting the physical associati on of the bright SMG s with 
the protocluster. In addition, lAlmaini et al. (2003) have 
reported strong clustering signals between SMGs and 
Chandra X-ray sources. These observational facts sug- 
gest that SMGs which host accreting massive black holes 
can tightly relate to large-scale structure at high-z, but 
the evidence for SMGs with heavily-obscured AGNs (i.e., 
protoquasars) forming in high-z protoclusters has not re- 
ported. 

Figured] shows the position of SSA22-AzTECl and the 
large-scale distributions of luminous objects towards/in 
the SSA22 protocluster. The high-density region of 
Lya emitters in SSA22 consists of at least three large- 
scale filamentary structures, and the highest density 
peak (i.e., the apparent protocluster core) actually cor- 
respon ds to the intersectio n of these large-scale fila- 
ments (|Matsuda et alJl2005h . SSA22-AzTECl is located 
close to the intersection of the filamentary structures, 
suggesting that this submillimeter-luminous protoquasar 
can trace the deepest central potential of the large-scale 
structure, possibly becoming a powerful AGN, just like 
high-z radio galaxies and quasars in protoc l usters (e.g. , 
McMahon et all 11991 iDunlop et all 11991 llsaak et al 



1994 1 : IIvisonlll995t lArchibald et al.l l2001t IStevens et al 



2003; IReuland et all 12001 IStevens et all 12010^ Al- 



though the accurate redshift of SSA22-AzTECl is yet 
to be determined, this may be the first observational 
example supporting the long-standing prediction that 
powerful AGNs preferentially form i n high-density envi- 
ronments such as protoclus ters (e.g.. iSiiacki et aU 120071 : 
iColberg fc Pi Matted [2001) . 
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Finally, we compare SSA22-AzTECl with other object 
with AGNs found in SSA22. Within the protocluster, 
at least over the Chandra coverage, no z = 3.1 AGNs 
that exceeds SSA22-AzTEC's AGN in X-ray luminos- 
ity have been reported. As shown in Figure |H1 how- 
ever, there are three X-ray-idcntificd AGNs at z = 3.1 
in the innermost part of the protocluster, whose rest- 
frame 0.5-8 keV lum i nosities (Lx ~ (2-3 ) x 10 44 ergs s~ 1 , 
iLehmer et alj|2009al : iGeach et al.ll2009() are comparable 
to that of SSA22-AzTECl. They have been identified in 
a Lya emitter and two Lya blobs in SSA22. 

An .SVi^ « 4-mJy source (SSA22-AzTEC6, 

iTamura et al.l [20091) could be associated with the X-ray 
luminous Lya emitter even though there is no submil- 
limeter interferometric identification. This implies that 
it might be another example of a protoquasar in the pro- 
tocluster. However, the effective photon index of the X- 
ray source (r e ff = 1.68^20) i s consistent with that of an 
unobscured AGN, and the host galaxy is clearly visible 
in the optical broad-band images unlike SSA22-AzTECl. 
These suggest that gas surrounding the AGN and stellar 
component is becoming to clear up. In contrast, the two 
Lya blobs within the high-density region are faint at 1.1- 
mirQ, which does not match the prediction for hosts of 
protoquasars that vigorous star-formation is associated 
with the early growth of massive black holes. The X-ray 
sources of the Lya blobs are unobscured (r c ff > 1), which 
may imply the absence of dense gas reservoir surround- 
ing the AGNs. Many studies predict that a luminous 
and unobscured phase of AGN evolution can follow the 
optically-thick phase at which the host is luminous in 
the FIR and the A GN is buried in dust clouds , as men- 
tioned above (e.g., [K^aSt^eOl][200l, [2006|). There- 
fore, we speculate that the AGNs in the Lya emitter and 
blobs are possibly at more advanced phase than SSA22- 
AzTECl's AGN. 

5. CONCLUSIONS AND FUTURE PROSPECTS 

Using the Submillimcter Array, we have identified an 
860-/im counterpart to an Si.i mm = 8.4-mJy source, 
SSA22-AzTECl, which has been previously discovered 
in a 1.1-mm single-dish survey towards the SSA22 pro- 
tocluster at z — 3.1. We then have investigated intrinsic 
nature of this source using multiwavelength data all the 
way from the radio to the X-ray. The derived properties 
of SSA22-AzTECl is summarized in Table [2] Our main 
findings and conclusions are as follows. 

1. We identified an 860-/im counterpart with a flux 
density of Sseofim = 12.2 ± 2.3 mJy and absolute posi- 
tional accuracy that is better than 0'.'3. The SMA coun- 
terpart is not resolved with the SMA 2" beam. The high 
brightness of this object is not due to strong gravitational 
lensing by foreground massive objects. (§ 13. ip 

2. The SED of SSA22-AzTECl is well constrained 
at wavelengths > 3 /zm, but drops out at a wavelength 
of 1 fim (J > 25.4 in AB, 2a, 1" aperture) or shorter. 
This steeply declining spectrum in NIR strongly suggests 
that SSA22-AzTECl can have a large obscuring column 
in front of the stellar component and/or an AGN. At 
the SMA position we see a significant X-ray counterpart 

23 A SCUBA 3.8 a detection has b een reported for one of the 
Lya blobs (LAB14, Gcach et al. 2005), but we find no significant 
emission (Si.imm < 2 mJy) in the better-quality AzTEC map. 



(?s20 counts), which has a hard X-ray spectrum (the pho- 
ton index, r eff = -0.34±g;||). (§ GO} 

3. It is important to determine the redshift of SSA22- 
AzTECl to see if it is located in the protocluster, but 
a conventional approach using rest-frame optical-to-NIR 
data does not provide a good constraint on redshift. 
On the other hand, photometric redshift estimates using 
MIR (24 /urn), submillimeter (860, 1100 /im), and radio 
(20 cm) data suggest the redshift of SSA22-AzTECl to 
be z = 3.19^0 35, consistent with the protocluster red- 
shift. (§E3J) 

4. We model the NIR-to-MIR SED of SSA22-AzTECl, 
and find that the NIR-to-MIR SED requires large extinc- 
tion (Ay ~ 3.4 mag) of starlight from a stellar compo- 
nent with M s tar ~ 1O 1O - 9 M assuming z « 3.1. (§ET5]) 

5. Modeling of the X-ray spectrum of SS A22- AzTEC 1 
suggests that it harbors a luminous AGN (Lx ~ 3 x 
10 44 ergs s _1 ) behind a large hydrogen column (Ah ~ 
1 x 10 24 cm~ 2 ) although it is in general hard to ac- 
curately determine Lx and Ah for AGNs with Ah < 
5 x 10 23 cm~ 2 . These observed properties of the AGN 
are very similar to those of the most heavily-obscured X- 
ray luminous SMGs that are likely candidates for proto- 
quasars. There are no robust candidates for protoquasars 
in high-z protoclusters other than SSA22-AzTECl, and 
it is thus possible that SSA22-AzTECl is the first exam- 
ple of a protoquasar growing at the core of a protocluster. 
(§§EHandE01) 

Spectroscopically determining the redshift for SSA22- 
AzTECl is particularly important to confirm the early 
growth of a massive black hole, as well as its host galaxy, 
in this high-density environment. Blind searches of CO 
molecular emission line(s) are necessary for spectroscop- 
ically identifying this object t hat is invisible in the op- 
tical and NIR. Zpectrometer (Harris et al.ll2007f) on the 
Green Bank Telescope, Z-Spec (jNavlor et al.l |200"3T ) on 
Caltech Submillimeter Observatory, EMI R on the IRAM 
30-m telescope (e.g.. IWeifi et all l2009afl . the Redshift 
Search Receiver (jErickson et al.ll2007l : lChung et alJ2009ft 
on the Large Millimeter Telescope, the 32-GHz wide- 
band spectrometer on the 45-m telescope at Nobeyama 
Radio Observatory, and the At acama Large Millime- 
ter/ submillimeter Array (ALMA, iWootten fc Thompson! 
120091 : ITguchi et al1l2009h will provide perhaps to achieve 
this. 
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Facilities: SMA, ASTE (AzTEC), VLA, Spitzer 
(IRAC, MIPS), Subaru (Suprime-Cam, MOIRCS), HST 
(ACS), GALEX (NUV, FUV), CXO (ASIS-I) 
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SMA identification of SSA22-AzTECl 



Table 1 

Multiwavelength counterparts to SSA22-AzTECl. 



Instrument 


Band 


Flux density 






Unit 


Ref 


VLA 


20 cm 


42 i 12 






MJy 
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Az 1 I Aa 1 Hi 


1 1 nn i , m 
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o q+0.8 a 
S -^-1.0 






mjy 


n 
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12.2 .: 2.7 






mjy 


1 


lvll.ro/ opiLZCT 


24 fim 


1 A 9 ^ 1S 

I^IZ IL lO 






MJy 


1 

1 


IRAC /Spitzer 


8.0 /im 


32.59 ± 2.27 








1 




5.8 fim 


22.24 ± 1.87 c 






fijy 


1 




4.5 fim 


14.46 ± 0.76° 






fijy 


1 




3.6 fim 


7.80 ± 0.50 c 






MJy 


1 


MOIRCS/Subaru 


2.14 fim (K s ) 


2.71 ± 0.63 c 






MJy 


1 




1.64 fim (H) 


2.13 ± 0.71 c 






MJy 


1 




1.26 /urn (J) 


< 1.12° 






MJy 


1 


S-Cam/Subaru 


905 nm (z') 


< 0.080 d 






MJy 


1 




771 nm («') 


< 0.042 d 






MJy 


1 




655 nm (H) 


< 0.032 d 






MJy 


1 




549 nm (V) 


< 0.035 d 






MJy 


1 




498 nm (iV,B497; 


i < 0.050 d 






MJy 


1 




448 nm (B) 


< 0.038 d 






MJy 


1 


ACIS-I/ Chandra 


0.5-8 keV 


2.71 c 


io- 


15 


erg cm -2 s - 1 


3 


0.5-2 keV 


9.27 c 


io- 


17 


erg cm~ 2 s~ 1 


3 




2-8 keV 


2.66 c 


10" 


15 


erg cm -2 s" 1 


3 



Refe rences. (1) This work; (2) ITamura et"aLl 120091) ; (3) ILehmer et"aLl 

(20093). 

a De-boosted flux density, i.e., flux density corrected for the flux bias due to confusion 

noise, with the 68% confidence interval. 

b 860-/jm flux density estimated from the visibility fitting. 

c Aperture-corrected flux density measured using a 3'.'0 aperture after the image is 
smoothed so that the PSF matches that of 8.0-/jm image. 
d la upper limit measured using a 2'.'0 aperture. 
e Flux integrated over the ACIS-I bands. 



Table 2 

The properties of SSA22-AzTECl. 



Coordinate (J2000) a z p hot b 
Right ascension Declination 


^FIR C 

(10 13 L o ) 


M star d 
(1O 1O M ) 


(ergs s ) 


22 h 17 m 32?42 +0° 17'44"01 3.19±g§| 


I Cj+0.4 


7 3+ 7 


« 3 x 10 44 



a SMA 860- fim position with an uncertainty of 0'.'33. 

b Photometric redshift estimate. The uncertainty is estimated from the 99% confi- 
dence interval. 

c Far-infrared luminosity. The uncertainty is estimated from the 99% confidence 
interval. 

d Stellar mass estimates derived from population synthesis models developed by 
IMarastofil 1120051 ). See text for the details. 

e The absorption corrected luminosity in rest-frame 0.5-8 keV. 
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Figure 1. {left) The SMA 860 /an continuum image towards SSA22-AzTECl (contours), overlaid on a Spitzer /IRAC false color image 
(8.0, 4.5, 3.6 /im for red, green, and blue, respectively). The SMA image is natural-weighted and CLEANed. The circle shows a field 
of view of the SMA (an FWHM of the primary beam, which is the point response function of the SMA's single dish elements) with a 
diameter of 34". The synthesized beam is indicated as an ellipse on the bottom- left corner. Significant emission 8" northward from the 
phase tracking center (cross) is evident. The peak intensity is 11.8 mjy beam -1 . The contours are drawn at (-4, -2, 2, 4, 6, 8)xcr, 
where lc = 1.4 mjy beam -1 . The negative levels are shown as dotted contours. In the inset of the left panel, we show the SMA 860-fim 
continuum image towards a test source J2218— 035. The systematic offset of the nominal peak position from the phase tracking center 
(cross) is (AR.A., ADecl.) = (O'.'OO ± 0'.'13, -0'.'31 ±ff.'X3). (right) The SMA 860 fim continuum image towards SSA22-AzTECl, overlaid 
on a Subaru/Suprime-Cam composite color image (blue, green, and red for B-, NB497-, and V-band, respectively). The green square 
denotes the position of the 'mini' Lyo blob LAB36, which is located ~ 10" away from SSA22-AzTECl. The circle shows a field of view of 
the SMA (34" in diameter). The contours are drawn as described above. 
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Figure 2. Visibility amplitudes versus projected baseline length for SSA22-AzTECl (solid circles). There is no evidence that SSA22- 
AzTECl is spatially resolved with the SMA beam. The visibility amplitudes are well fitted to a constant value 12.2 ± 2.7 mjy (solid line) 
as a function of projected baseline length. The 1<t confidence interval is shown as a region sandwiched between two dotted lines. The 
phase center is shift to the peak position of the 860-/xm counterpart to SSA22-AzTECl. We also show visibility amplitudes expected for 
axisymmctric 2-dimensional Gaussians with FWHM of 0"5, l'.'O, and 2'.'0. 
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Figure 3. Postage stamp images of SSA22-AzTECl at radio to near-ultraviolet wavelengths. The red circles mark the SMA position. 
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Figure 4. The Spitzer /IR AC Ss.g/ S3sSg.o/ S4.5 color-color diagram for SMGs. The rilled circles in blue, green, and red are SMGs with 
millimeter CO or Spitzer/IRS MIR spectroscopy at z < 1.5, 1.5 < z < 3.0, and z > 3.0, respectively. The big red star is SSA22-AzTECl, 
and indeed it is found near t he transition regi on between 2 < 3 and z > 3 SMGs. The smaller, black stars represent infrared-luminous 
quasars from the FLS survey {Lacy et al. 2004), and it seems th e z > 3 SMGs ind eed have colors that overlap with the infrared power-law 
AGNs (thick magenta line). The color region on the plot that Yun ct al. (200||) proposed for the identification of SMG counterparts is 
outlined by dash-dotted line. The region for infrared-luminous AGNs proposed by[Lacv et al. (2004) is outlined by the long dashed line. 




Redshift Wavelength (/im) 

Figure 5. Radio-to-MIR photometric redshift of SSA22-AzTECl. (a) The top panel shows the reduced-^ 2 versus redshift. Template 
SEDs arc from Michalowski ct al. (2010). The grey curves show the reduced-x 2 for 35 SEDs fit to 24-/im-dctcctcd SMGs with spectroscopic 
redshift. The thick solid curve shows that for an SED averaged over the 35 SEDs, and the horizontal dashed line shows the 99% confidence 
interval. The yellow bar represents the redshift range of the protocluster (2 = 3.06-3.12). The bulk of the curves have their local minima 
at the protocluster redshift, and hence we adopt z = 3.1 in the rest of the an alyses. The bottom panel shows the infrared luminosities that 
allow the mi nimum y 2 value at each redshift. The details are described in § 13.4.21 (b) Th e observed flux densiti es (filled circles) and the 
SED (curve, Michalowski ct al. 2010) averaged over 35 SMGs with spectroscopic redshift (Chapman ct al. 20051). The SED is shifted at 
the best-fit redshift (2 = 3.19) and is scaled to the best-fit FIR luminosity (logLpiR = 13.3). Only the data at A t, s > 24 fim are used in 
the x 2- fit (open circles). The thin curve indicates the best-fit SED model (Maraston 2005) for the stellar component of SSA22-AzTECl. 
The rapid dimming and break at A t, s < 8 ^m is likely due to obscured stellar (a nd p ossibly AGN) emission that suffers from extremely 
heavy attenuation (Ay fs: 3.4) of gas and dust in SSA22-AzTECl. See the text (§ [375J for more details. 
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Figure 6. The position of SSA22-AzTECl and the large-scale distributions of lu minous objects towa rds/in the SSA22 protocluster; 
AzTEC sources (fill ed circles, whose size s are proportional to 1.1-mm flux densities; Tamura ct al. 2009), robust candidates for z = 3.1 
Lyct blobs (crosses; Matsuda ct al. 2004), and the most luminous z = 3.1 AGNs in SSA22 with rest 0.5—8 keV luminosities of Lx ~ 
(2-3)xl0 44 ergs s" 1 (squares; iLehmer et 'al. 2009a; Gca ch et al.ll2009T) . The white square denotes the heavily-obs cured AGN of SSA22- 
AzTECl. The gray-scales show the projected number density of Lya emitters at z = 3.1 (Havashino ct al. 2004) and the contours are 
drawn at (1, 1.5, 2.0, 2.5) times the mean density of the Lya emitters, which outlines the z = 3.1 protocluster. The dashed line shows the 
field-of-view of the Chandra Deep Protocluster Survey. The bar at bottom-right corner indicates a comoving scale of 10 Mpc at redshift of 
z = 3.1. 



